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A review of research on the effects of gravity on combustion processes is presented, with an emphasis 
on a discussion of the ways in which reduced-gravity experiments and modeling has led to new under- 
standing. Comparison of time scales shows that the removal of buoyancy-induced convection leads to 
manifestations of other transport mechanisms, notably radiative heat transfer and diffusional processes 
such as Lewis number effects. Examples from premixed-gas combustion, non-prenuxed gas-jet flames, 
droplet combustion, flame spread over solid and liquid fuels, and other fields are presented. Promising 
directions for new research are outlined, the most important of which is suggested to be radiative reab- 
sorption effects in weakly burning flames. 


Introduction 

Gravity influences many combustion processes, 
particularly due to buoyant convection which affects 
transport of thermal energy and reactants to and 
from the chemical reaction zones. Recently, many 
experimental and theoretical studies of combustion 
at microgravity (jug) conditions have been con- 
ducted. These studies are motivated by the need to 
assess fire hazards in spacecraft and to enable better 
understanding of combustion processes at earth 
p*avity (lg) through the elimination of buoyancy 

This paper discusses how new understanding of 
combustion processes has been obtained through fig 
research, rather than providing a comprehensive 
review of this rapidly changing field. First, com- 
parisons of time scales for various chemical and 
transport processes in flames, including buoyancy- 
induced transnort, are given. Next, examples of un- 
expected results and new understandings obtained 
through fug research are discussed. These findings 
are then summarized and future research directions 
are suggested. 

Comparison of Time Scales for Premixed-Gas 
Combustion 

To determine the conditions in which gravity can 
affect flames, the estimated time scales for chemical 
reaction (t chem ), inviscid buoyant convection (f Snv ), 
viscous buoyant convection (f^), conductive heat 
loss to walls (tcond), and radiant heat loss (t rad ) are 
compared. Premixed laminar flames are considered 
first because of their simplicity. Subsequent sections 
introduce time scales for other flames. 


The chemical time scale is (see Nomenclature) 
fchem ~ <*/S L where S = a/S L , thus t cliem - a/St- The 
convective transport time scale is dJU, where d is a 
characteristic flow length scale, U ** ( gd(Ap/p)) v 2 is 
the buoyant convection velocity, and Ap the density 
change across the flame. Because Aptp ** 1 for 
flames, t inv ** d/(gd) l/2 = ( dJg ) v 2 . For inviscid flow, 
d is determined by the apparatus dimensions, for 
example the burner or tube diameter. For viscous 
flow, d cannot be specified independently; instead d 
v/U, thus U ** (gv) 1/3 and t^ s ** dJU (v/U)/U *■ 
(v/g 2 ) l/3 the conduction time scale f a)nd is the flame 
front temperature (7» divided by the rate of tem- 
perature decrease due to conductive loss; thus f tond 
- T { /(dT/dt) ~ T ( /(pCMTf - TJ); thus - 
d 2 /! 6a. Similarly, for optically thin radiation, f rild 
T[/(A/pC P ) - f y/(y - l)|(P/4 a a? (T* - T*)). 

Two sets of time scales are shown in Table 1, one 
for near-stoichiometric hydrocarbon-air flames and 
one for near-limit flames, both at P = 1 atm. For 
near-stoichiometric flames S L ** 0.40 m/s, T ( ** 2200 
K, a ** v 1.5 X 10“ 4 m 2 /s and a ? ** 0.56 m“ l . For 
near-limit flames, S L 0.02 m/s, T f 1500 K, a *- 
v ^ 1.0 X 10 -4 m 2 /s and a P *** 0.83 m _l . For both 
cases g ■» 9.8 m/s 2 , y *** 1.35, T* 300 K and d 

0. 05. m (a typical apparatus dimension.) 

Several observations can be made on the basis of 
these simple estimates: 

1. Buoyant convection is unimportant for near-stoi- 
chiometric flames because f Nls > t chem and f hlv > 

^chetn • 

2. Buoyant convection strongly influences near-limit 
flames at lg because < f chen) and t m , < t chem . 

3. Radiation effects are unimportant at lg compared 
to buoyant convection because f vis f rjd and f inv 
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TABLE 1 


Estimates ol time scales tor stoichiometric and near- 
limit hydrocarbon-air Haines at l atm pressure 


Time Scale 

Stoichiometric 

Flame 

Near-Limit 

Flame 

Chemistgry (7, | ((M J 

0.00094 s 

0.25 s 

Buoyant, inviscid (f im ) 

0.071 s 

0.071 s 

Buoyant, viscous (f v1 j 
Conduction to tube wall 

0 012 s 

0.010 s 

i) 

0.95 s 

1.4 s 

Radiation (f r . ld ) 

0.13 s 

0.41 s 


4. Radiation effects dominate near-limit flames be- 
cause f rac i ** t L . | K . ni , but these effects are only ob- 
servable at reduced gravity' because of observa- 
tion 3. 

5. The apparatus size (d) must be larger than about 
0.03 m to observe radiation-induced extinction; 
otherwise, conduction losses exceed radiative 
losses (f Hjnd < f rad ). 

6. Many radiative loss effects can be studied in 
drop-towers (test duration 2—10 s), because these 
times are typically larger than t rad . 

7. Because t im — g l/2 and t^ s — g 1/3 , aircraft-based 
Mg experiments at g 10 ' 2 g 0 may not provide 
sufficiently reduced buoyancy to observe radia- 
tive effects. 

8. Because f vi5 - v 1/3 - p- and t rad - p/A - 
P'/P [ ~ P°, t vis /f rad ~~ P~ 1/3 Thus, t rad is indepen- 
dent of P , but at higher P, buoyancy effects in- 
terfere more strongly with radiative effects. 

9. A Reynolds number Re d » Ud/v for buoyant flow 
is estimated as (gdy^d/v 2 = Cr 1 / 2 , where Gr d = 
gcP/v 2 is a Grashof number. For Re d > 10 3 , thus 
Gr d >: 10 6 , buoyant flow at Ig is generally tur- 
bulent, thus, it difficult to obtain steady laminar 
flames in large systems at lg. 

The implications of these observations are dis- 
cussed in the following sections. Essentially, any 
combustion process where f L . hein or t rad exceeds f inv 
or t vis may be affected by gravity and is worthy of/*g 
investigation. 


Premixed Gas Flames 


since limits are different for upward, downward, ancl 
horizontal propagation [5]. 

Practically all flammability' limit studies show that 
burning velocity at the flammability limit IS, JtlM ) is 
nonzero. Giovangigli and Smooke [6] have shown 
that there is no purely chemical flammability limit 
criterion for planar unstretched flames; without 
losses, S L decreases asymptotically to zero as dilution 
increases. Consequently, loss mechanisms such as 
those discussed below are needed to explain limit 
mechanisms. The resulting predictions of S Uim in- 
dicate that S Ultn usually depends only vveaklv on 
chemical reaction rate parameters. Thus, limit 
mechanisms may be inferred bv comparing pre- 
dicted and measured S Uhn without detailed chemi- 
cal knowledge. The mixture composition at the limit 
affects S Uim only weakly through 7>: thus, compar- 
ing predicted and measured limit compositions is not 
especially enlightening; comparisons of S Llim values 
is much more useful. Consequently this discussion 
emphasizes comparisons of predicted and measured 
values of S Uim . 

For upward propagation, Levy [7] showed that the 
flame rise speed at the limit 0.33(gd) 1/2 ) is iden- 
tical to that of an inviscid hot gas bubble. This re- 
lation was later verified for a wide range of tube di- 
ameters and mixtures [8]. Buckmasterand Mikolaitis 
[9] showed how this minimum rise speed causes hy- 
drodynamic strain at the flame tip, which causes ex- 
tinguishment for sufficiently low S L . The predicted 
burning velocity of the limit mixture (S L j im ) is, after 
temperature-averaging transport properties: 


^L.lim 2.8 exp 




The form of equation 1, S L!jm - ( gcr/d ) u 4 , can be 
obtained by setting t inv = t chenv 
For downward propagation in tubes, centrifuge 
experiments [10] indicate that S Uim - g l/3 , inde- 
pendent of Le , which is reasonable since downward- 
propagating near-limit flames are nearly flat and 
unstrained. Experiments [11] and numerical simu- 
lations [12] suggest extinction results from sinking 
regions of cooling burned gas near the walls over- 
taking the flame and blocking it from fresh reactants, 
although the g 1/3 scaling was not tested in Refs. 
[11,12]. The g^ 3 scaling can be obtained by setting 
^chem ^is> thus 


Flammability Limits 

The previous section showed that gravity effects 
are significant only for mixtures with low S L , imply- 
ing mixtures highly diluted with excess fuel, oxidant, 
or inert gas, but sufficient dilution causes flamma- 
bility limits. Thus, gravity effects have significant in- 
fluences on near-limit behavior, which is expected 


Sulim - (ga) l/3 (2) 

Experiments [8] employing varying diluent gases 
and pressures confirm the a 173 scaling and lack of 
dependence on tube size, which, together with the 
g 173 scaling mentioned above, supports the proposed 
mechanism. 

Both upward- and downward-limit mechanisms 
indicate that as g — > 0, S L ,Um 0, implying arbitrarily 
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Fic. 1. Predicted values of burning velocity and peak 
flame temperature in CH 4 - (0.21 0 2 4* 0.49 N 2 + 0.30 
C0 2 ) mixtures under adiabatic conditions, with optically 
thin radiative losses and including reabsorption effects 
(26]. 


weak mixtures could bum very slowly. However, 
conductive or radiative losses prevent arbitrarily 
weak mixtures from burning even at g = 0. Theories 
that relate flammability limit to heat losses [13-15] 
predict a minimum Tf below which flame propaga- 
tion cannot occur because chemical reaction rates 
are much stronger functions of temperature than 
heat loss rates (exponential vs. algebraic). Conse- 
quently, because dilution decreases Tf, dilution in- 
creases the impact of losses, leading to flammability 
limits. For conductive losses, setting f chem = f colld 
leads to 


f^iim = S Uini d/a = constant (3) 

with experiments [8,16] and computations [17] in- 
dicating Pe hm 40. For radiative losses, setting f chem 
= t rad leads to [18,19] 


P*C p 


1 . 2 PAX< 


(4) 


For lean-limit CH 4 -air mixtures at 1 atm, equation 
4 yields S Llim ** 0.023 m/s, similar to detailed nu- 
merical model predictions [20,21]. Such small S L liin 
are not observed at Ig because of buoyant convec- 
tion (f jnv < f ra d and f vis < f rad ); equations 1 and 2 
yield j im *» 0.033 and 0.078 m/s for upward and 
downward propagation, respectively. At /ig, however, 
predictions of equation 4 compare favorably to ex- 
periments in large combustion vessels [19,22] using 
varying pressures, fuels, and inert gases. Also, similar 
results were obtained for CH 4 -air mixtures at 1 atm 
in tubes with d = 0.05 m [23], suggesting these lim- 
its are apparatus-independent. Thus, radiative losses 
may cause flammability limits when extrinsic losses 
(conduction, buoyant convection, etc.) are elimi- 
nated. In this instance, /ig experiments enabled ob- 
servation of phenomena not observable at lg. 


These radiative effec ts apply only for optically thin 
gases (i.e., no reabsorption of emitted radiation), 
which are inappropriate for large systems, high pres- 
sures, or mixtures with strongly absorbing material. 
With this motivation, fig experiments [24 [were con- 
ducted using lean CH 4 -air mixtures seeded with SiC 
particles. Because solids eimt/absorb as black- or 
grav-bodies, whereas gases radiate in narrow spectral 
bands, particle-seeded gases emit/absorb more ra- 
diation than particle-free gases. Measurements of 
propagation rates, pressures and postflame thermal 
decays showed that, consistent with theoretical pre- 
dictions [25], at low particle loadings the particles 
increase radiative loss (optically thin conditions), 
whereas at higher loadings reabsorption of emitted 
radiation becomes significant, which decreases net 
radiative loss and augments conductive heat trans- 
port. 

Even for gases, computations [26] using detailed 
statistical narrow-band radiation models show that 
flammability limits are extended considerably with 
reabsorption (Fig. 1). With gases, however, two 
mechanisms lead to flammability' limits even with 
reabsorption. One is the difference in composition 
between reactants and products; if H 2 0 or other ra- 
diatively active combustion products are absent from 
the reactants, radiation from these species that is 
emitted upstream cannot be reabsorbed by the re- 
actants. The second mechanism is that emission 
spectra are broader at T f than T x ; thus, some radia- 
tion emitted near the flame front cannot be absorbed 
by the reactants. Upstream loss occurs via both 
mechanisms, leading to extinction of w>eak mixtures. 
These results suggest that fundamental (domain- 
and gravity-independent) flammability limits due to 
radiative losses may exist at /ig, but these limits are 
strongly dependent on emission/absorption spectra 
of reactant and product gases and their temperature 
dependence ana cannot be predicted using gray-gas 
or optically-thin model parameters. 


Stretched Flames 

Premixed gas flames generally are not flat and 
steady nor do they propagate into quiescent flows. 
Consequently, flames are subject to “flame stretch," 
2T = (1 IA)(dA/dt) [18], which affects S L and extinc- 
tion conditions [18,27]. At lg, buoyancy imposes 
flame stretch comparable to t " v ! or t ~ s l . At /ig, weak 
flame stretch effects that are insignificant at lg may 
dominate. One example is expanding spherical 
flames for which 


1 — 
A dt 


1 d . , 2 dr ( 

- — s -r (4niy) = - — 

4 nrj dt 1 r dt 


(5) 


For Le < 1, positive stretch increases T ( because the 
increased chemical enthalpy diffusion to the flame 
front in the form of scarce reactant exceeds the in- 
creased thermal enthalpy loss. Because heat release 
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Fig. 2. Characteristics of self-extinguishing flames in 
CH 4 -air mixtures at 1 atm for various mole percent CH 4 
and spark ignition energies [112]. The "x” symbols denote 
extinction. 



Fig. 3. Measured and predicted extinction strain rates 
for strained premixed CH 4 -air flames at fig [33] showing 
dual-limit behavior, that is, residence-time limited extinc- 
tion at high strain rates (upper branch, "strong flames”) and 
radiative loss extinction at low strain rates (lower branch, 
“weak flames”). 


reactions have high activation energies, small T( 
changes cause large changes in reaction rate and 
thus S L . An evolution equation for nonadiabatic ex- 
panding spherical flames is given by [28]: 

S +S21nS2 = f- <? (6) 

where S - (dr { /dt)/(S L (pJp ( )), R * r ( !(fi31(Le, e)). 
I(Le, e ) is a scaling function (/ > 0 for Le < 1 and 
1 < 0 for Le > 1) and Q » {fiA( T T )fi)/{X(T f - TJl 
The terms in equation 6 represent unsteadiness, 
heat release, curvature -induced stretch, and heat 
loss, respectively. For steady planar flames, equation 
6 becomes S 2 In S z = - Q , which exhibits a maxi- 
mum Q = 1/e = 0.3678 . . . at S « e~ 1/2 f which 
corresponds to S L ^ from equation 4. For Le < 1, 


the curvature effect (2 S/R) opposes heat loss (£), 
allowing mixtures that are nonflammable as plane 
flames (£> > 1/e) to exhibit expanding spherical 
flames until r, grows too large and thus the curvature 
benefit too small. For mixtures just outside the limit, 
the extinction radius may be very' large. Such behav- 
ior, termed self-extinguishing flames (SEFs), is ob- 
served experimentally [19,22] (Fig. 2) a t fig j n near- 
limit mixtures with Le slightly less than unity. 
( Mixtures with lower Le exhibit diffusive- thermal in- 
stabilities or flame balls discussed below.) Equation 
6 also predicts, consistent with experimental obser- 
vations [29], that SEFs cannot occur for Le > 1 
(thus, R < 0) because both cur/ature and heat loss 
weaken the flame. 

Two experimental observations not predicted by 
equation 6 are that narrow mixture ranges exhibit 
both SEFs and normal flames and that the energy 
release before extinguishment can be orders of mag- 
nitude greater than the ignition energy. Such behav- 
ior is predicted by computations [30] not subject to 
scaling limitations of activation energy asvmptotics 
used to derive equation 6. These calculations also 
show that for small initial r f , all mixtures exhibit ex- 
tinguishment, corresponding to nonignition behav- 
ior [31]. Thus, in mixtures exhibiting SEFs, flames 
extinguish at large curvature (small r f ) due to large 
27 and at small curvature due to radiative losses. This 
dual-limit behavior is also exhibited by many other 
types of flames described later. 

Flames in hydrodynamic strain induced by coun- 
terflowing round jets are frequently employed to 
model turbulence-induced flame stretch effects. At 
steady state, the flame resides at the axial location 
iy) where the axial velocity (U y ) equals S L for the 
given 27 = dU,/dy. As L increases, U y increases, thus 
the flame moves toward the stagnation plane 
(smaller tj) and the burned gas volume (thus radia- 
tive loss) decreases. As with curvature-induced 
stretch, for Le less than/greater than unity, moderate 
hydrodynamic strain increases/decreases S L , but for 
all Le , large strain extinguishes the flame [32], Con- 
sequently, fig experiments [33] in low-Le mixtures 
(Fig. 3) reveal extinction behavior analogous to 
spherical flames. For large 27, the short residence 
time (~2T _1 ) causes extinguishment (27' 1 - f chein ) 
along the "normal flame” branch, analogous to non- 
ignition behavior of spherical flames. In contrast, for 
low Z, the residence time and burned gas volume 
are large; thus, radiative loss is significant (f rad 
*chem)> so radiative loss extinguishes the flame along 
the “weak flame” branch, analogous to SEFs. The 
optimal Z (**13 s* 1 ) producing the minimum flam- 
mable fuel concentration corresponds to 27* 1 = 
0.08 s, which is less than t yi5 or f inv ; thus, the C- 
shaped response and the entire weak-flame branch 
cannot be observed at lg. The optimal 27 is nearly 
the same for model and experiment, indicating that 
loss rates are modeled well, but the computed limit 
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Fig. 4. Schematic diagram of a flame ball, illustrated for 
the case of fuel-limited combustion at the reaction zone. 
The oxygen profile is similar to the fuel profile except its 
concentration is nonzero in the interior of the ball. The 
combustion product profile is identical to the temperature 
profile except for a scale factor. 



Fig. 5. Predicted effect of heat loss on flame-ball radius 
and stability properties [42] showing radially unstable 
(small) flame-ball solution, radially stable (large) flame-ball 
solution, and three-dimensional instability for large flame 
balls. 


composition is leaner than the experimental limit, 
suggesting that the chemical mechanism used is in- 
accurate for weak mixtures. Due to the radiant loss 
decrease at moderate 27, the flammability limit ex- 
tension also occurs for Le > 1, though, for suffi- 
ciently high Le, no C-shaped response or flamma- 
bility limit extension occurs [34]. (For spherically 
expanding flames, no flammability extension occurs 
for Le > 1 because in this case there is no mecha- 
nism to reduce radiative loss by flame stretch.) 

The combination of nonmonotonic response to Z 
plus the reduced radiative loss at larger 27 causes 
several new extinction branches depending on 


t v | lrm , t fM j, 27' l , and Le [34,35]. It is uncertain 
whether these branches are physically observable 
because they have not been identified experimen- 
tally and stability’ analyses have not been performed. 


Flame Balls 

Over 50 years ago, Zeldovich [36] showed that the 
steady mass, energy, and species conservation equa- 
tions admit solutions corresponding to stationary 
spherical flames, characterized by a flame radius (ly ). 
Fuel and oxygen diffuse from the ambient mixture 
inward to the reaction zone while heat and combus- 
tion products diffuse outward (Fig. 4). Mass conser- 
vation requires that the fluid velocity be zero every- 
where. The temperature and species mass fraction 
profiles have the form c j + cJr, where c x and c z are 
constants. Corresponding solutions in planar and cy- 
lindrical geometry cannot exist because the solution 
forms Ci + c 2 r and c x + c 2 In(r), respectively, are 
unbounded as r — » oo. Zeldovich [36] and others 
[37,38] also showed that flame-ball solutions are un- 
stable and thus probably not physically observable, 
although these solutions are related to flame ignition 
[37]. 

Forty years after Zeldovich [36], apparently stable 
flame balls were accidentally discovered in drop- 
tower experiments using H 2 -air mixtures [39] and 
aircraft-based fig experiments using various low -Le 
mixtures [40]. The/^g environment facilitated spher- 
ical symmetry and prevented buoyancy-induced ex- 
tinction. For mixtures sufficiently far from flamma- 
bility limits, expanding spherical fronts composed of 
many individual cells were observed, whereas for 
more dilute mixtures, cells that formed initially did 
not split and instead closed up on themselves to form 
flame balls. (For still more dilute mixtures all flames 
eventually extinguished.) It was inferred that flame 
balls can occur in all near-limit low -Le mixtures; 
however, the short duration of drop-tower experi- 
ments and substantial g fluctuations in the aircraft- 
based fig experiments precluded definite conclu- 
sions. Recent space shuttle experiments [41] 
confirmed that flame balls can exist for at least 500 
seconds (the entire experiment duration.) 

Zeldovich [36] noted that radiative losses might 
stabilize flame balls; consequently, after their exper- 
imental observation, radiative loss effects on flame 
balls were analyzed [42]. For moderate loss, two so- 
lutions branches are predicted (Fig. 5), a strongly 
nonadiabatic large-raaius branch and a nearly-adia- 
batic small- radius branch. For sufficiendy strong 
losses, no solutions exist, indicating extinction limits. 
Stability analyses [42] predict that all small flames 
are unstable to radial disturbances, and large flames 
with weak loss (far from flammability limits) are un- 
stable to three-dimensional disturbances. Close to 
extinction limits, the large- radius branch is stable to 
both disturbances. These predictions are consistent 
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Fig. 6 Comparison of computed flame-ball radii as a 
function of Hi mole fraction in Hi-air mixtures for three 
different H 1 -O 1 chemical mechanisms, along with prelim- 
inary results from the STS -83 and STS-94 space experi- 
ments [45]. 


1 6 | 1 1 1 1 1 r 



2 4 6 8 10 12 14 16 

Mole % H 

z 

Fig. 7. Computed flame-ball radii as a function of the 
H 2 mole fraction for steady flame balls in H 2 - 02 -C 0 2 
mixtures with H 2 :0 2 = 1:2, for optically thin C0 2 radiation 
and with COi radiation artificially suppressed (optically 
thick limit for C0 2 radiation). Preliminary experimental re- 
sults from the STS-94 mission are also shown (filled circles) 
[45J. 

with the observed splitting cellular flames away from 
limits and stable balls close to limits. For Le close to 
or larger than unity, all flame balls are unstable for 
any loss magnitude [43], explaining why they are 
never observed in (for example) CH^-air mixtures 
(Le 0.9) or C^Hq-air mixtures (Le ** 1.7). 

Numerical predictions of non-adiabatic flame balls 
employing detailed chemistry, diffusion, and radia- 
tion models [44,45] are qualitatively consistent with 
these experimental and theoretical results. Still, 
quantitative agreement has been elusive (Fig. 6) for 
at least two reasons. First, flame-ball properties are 
very sensitive to the three-body recombination step 


H + O? + HiO -4 HO-> + HiO [45] whose rate 
varies widely between different published H->-Oi re- 
action mechanisms. The second reason is reabsorp- 
tion of emitted radiation in mixtures diluted with 
radiatively active C0 2 or SF H . An upper bound on 
self-absorption of diluent radiation ( a H -4 ^c) is as- 
sessed by neglecting diluent radiation entirely be- 
cause as rz P 2o, radiative loss from the diluent van- 
ishes and no additional heat transport occurs due to 
radiation. Agreement between predicted and mea- 
sured flame radii is much better in this case (Fig. 7), 
strongly suggesting radiation modeling including 
reabsorption is needed for accurate predictions in 
these cases. 

A key difference between propagating flames and 
flame balls is that propagating flames have convec- 
tive-diffusive zones where temperature and concen- 
tration approach their ambient values in proportion 
to e~ r/t \ whereas flame balls have purely diffusive 
zones where the approach is proportional to 1 Ir. 
Plane flames respond on short time scales t c i u . 1T1 — 
6 2 /a, whereas the gradual 1/r flame-ball profiles pro- 
duce properties dominated by the far- field length 
scale firf, and thus diffusion time scales (/fry ) 2 /a [42], 
typically 100 s. These scales are relevant to stability 
and extinction limits because they affect the times 
for radiant combustion products to diffuse to the far- 
field and indicate large volumes of gas ( ~~(Pr } ) 
where radiative loss affects flame balls. Such large 
scales are confirmed by space experiments [41] and 
numerical simulations [44,45]. Droplet and candle 
flames (discussed later), which have quasi-spherical, 
diffusion-dominated far-fields, exhibit analogous be- 
havior. 


Caseous Non-Premixed Flames 

Stretched Flames 

Non-premixed flames, where fuel and oxidant are 
separated before combustion, are affected by stretch 
differently than are premixed flames. The most sig- 
nificant difference is that the flame position is de- 
termined by the location of stoichiometric mixture 
fraction, dictated by mixing considerations, rather 
than being determined by balances between S L and 
U as in premixed flames. Consequently, non-pre- 
mixed flames have considerably less freedom of 
movement. Also, premixed flames have characteris- 
tic thicknesses 6 ~ a/S L unrelated to the flow envi- 
ronment, whereas non-premixed flames have only 
the diffusion length scale 6 — (a/E) m . With fixed 
flame location and 5 increasing monotonically with 
decreasing 2T, non-premixed flames with radiative 
loss exhibit only simple C-shaped responses to strain 
(Fig. 8) [46], with a short residence time extinction 
branch (f c h em > £~ l ) and a radiative loss extinction 
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Fig. S. Measured and predicted extinction strain rates 
for strained non-premixed X r diluted CH 4 versus air coun- 
terflow flames at fig [46] showing dual-limit response anal- 
ogous to preniLxed flames (Fig. 3). 


TABLE 2 

Predicted scalings of flame heights (L t ) and residence 
times (/^ t ) for non-premixed round-jet and slot-jet flames 
under momentum -dominated and buoyancy- 
dominated conditions 


Geometry 

Flow 

Mechanism 

4 

tje, 

Round-jet 

Momentum 

U,4f/D 

<iyo 

Round-jet 

Buoyant 

U^f/D 

(U^ygD)^- 

Slot-jet 

Momentum 


df/D 

Slot-jet 

Buoyant 

{UffiJD* g) w 



branch (f rac j < L~ *) rather than the compli- 
cated responses found for premixed flames [34], The 
only significant difference in flame structures near 
the two limits is S [46], a situation quite unlike pre- 
mixed flames. 

For the radiative extinction branch f c j, ein is still a 
factor because f rat j ** 1 results in order unity de- 

creases in flame temperature, thus causing exponen- 
tially large decreases in f c j lem . Even conditions far 
from extinction at lg, therefore, may exhibit radia- 
tive extinction at fig due to much larger residence 
times. This mechanism also applies to radiative ex- 
tinction of other types of non-premixed flames dis- 
cussed later. It is also somewhat analogous to the 
lower branch of strained premixed flames (Fig. 3) 
and the large-radius branch of flame balls (Fig. 5). 

The experiments shown in Fig. 8 suggests that no 
flames exist below some value of L, whereas the 
model predicts flames at arbitrarily low 27. Similar 
behavior was seen for premixed flames (Fig. 3). This 
suggests an additional loss mechanism not consid- 
ered by the model, probably axial conductive heat 


losses to the jets or radial conductive loss to inert 
gases surrounding the reactant streams. This would 
induce f t J M | ** 2.9 s ~ 1 if cl is the jet spacing (25 mm), 
or 7. 1 s ~ 1 if d is the jet diameter ( 16 mm ), Either of 
these are roughly consistent with the minimum 27 in 
Fig. 8. Thus, apparatuses large enough to studv 
flames at lg without substantial conductive loss, 
w’here maximum length scales are about ( of vis ) 12 , are 
insufficient for the weaker flames attainable at fig, 
where maximum length scales are about 


Laminar Gas-Jet Flames 


A fuel jet issuing into an oxidizing environment is 
one of the simplest types of flames. Jost [47] and 
Roper [48] estimated the flame height (Lf) and res- 
idence time from jet exit to flame tip (f jt . t ) by deter- 
mining the height ( tj) where the transverse diffusion 
time ( d(y) 2 /D , where d(y) is the stream tube width) 
equals the convection time (U(y)/y, where U(y) is 
the axial velocity). When buoyancy and viscosity ef- 
fects are negligible (momentum-controlled jets), 
U(y) is constant and equals the jet exit velocity (U,X 
whereas when buoyancy effects dominate, U(tj) ~ 
(gy) 1/2 . In either case, mass conservation requires 
that d(i/) 2 £/(i/) = d 2 (7 0 = constant for round jets or 
d(y)U(y) = d 0 U 0 = constant for slot jets. The re- 
sulting estimated scalings for L f and f jet are given in 
Table 2. Transition from buoyancy-controlled to mo- 
mentum-controlled conditions occurs when the time 
scale for the former exceeds the latter, which cor- 
responds to U Q > gdfyD for either round jets or slot 
jets. The scalings for momentum-dominated flames 
presume constant U, which is reasonable for coflow- 
ing Burke-Schumann flames, but for nonbuoyant jet 
flames without coflow, the jet spreads and deceler- 
ates. For this situation [49] 


Li 


UJl 1 

D 2Sc l/ * 



= U Jl. = l n (_ 

l U(y) D \l 



(7) 


Because Sc 1, the scalings of L r and fj et are similar 
with or without viscosity. Fig. 9 shows measurements 
of Lf for CH 4 flames [50]. Note that, as the scalings 
predict, L(/d 0 — Re 55 U 0 dJ v at both lg and fig, and 
only small differences exist between lg and/ig flame 
lengths. 

All fig studies show larger flame widths (w) at fig 
than lg due to lower U and longer [50]. Also, w 
is larger at fig because the temperatures are lower 
(see below) and D ~ T 1 Because w depends on 
whether C7 is accelerating (buoyant jets), constant 
(nonbuoyant Burke-Schumann flames), or deceler- 
ating (nonbuoyant jets), w is more difficult to predict 
than is Lf [2]. The difference between lg and fig 
widths decreases as Re (thus, U 0 ) increases (Fig. 10). 
The nonbuoyant widths increase slightly with Re, 
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Fic. 9. Measured flame lengths, normalized by jet di- 
ameter, as a function of the jet Reynolds number for non- 
p re mixed CH 4 -air jet flames at lg and fig. Data is taken 
from a variety of sources and compiled in Ref. [50 J. The 
data show a nearly linear relationship between flame length 
and Reynolds number, with generally longer flame lengths 
at/*g, due to the differences between residence times un- 
der buoyancy-driven versus momentum -driven residence 
times. 
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Fig. 10. Measured flame widths, normalized by jet di- 
ameter, as a function of the jet Reynolds number for non- 
premixed CH 4 -air jet flames at lg and fig. Data is taken 
from a variety of sources and compiled in Ref. [50]. The 
data show larger flame widths at fig due to the differences 
between accelerating flow at lg versus decelerating flow at 
fig. The data also show that, consistent with theoretical pre- 
dictions, the width is nearly independent of Reynolds num- 
ber for nonbuoyant conditions, except at low Re where 
boundary- layer approximations are invalid. 

whereas all aforementioned models predict self- 
similar flame shapes with no effect of Re. This may 
be due to axial diffusion, not considered in these 
models, which increases mixing over that with radial 
diffusion alone. This suggestion is supported by the 
results shown in Fig. 10 which show that w/d n is 
lowest for lowest Re (^SO), where axial diffusion is 



Fig. 11. Measured flame heights for non-premixed 
C 3 H s -air jet flames at lg and/ig [5l| showing transition to 
turbulence. Note lg flame lengths are shorter than fig 
flame lengths, even at very high Reynolds numbers. 


most significant, but asymptotes to fixed values at 
high Re. 

Compared to yellow lg flames, fig gas-jet flames 
are more red [49,50], indicating lower blackbodv 
soot temperatures and presumably lower maximum 
flame temperatures. This occurs because is larger 
at Mg and thus radiative loss effects are 

greater. Drop-tower [51] and space [52] experiments 
indicate surprisingly large and consistent radiative 
loss fractions (0.454).60) at jug compared to 0.07- 
0.09 at lg, for various fuels, pressures, 0 2 mole frac- 
tions, and flow rates. Thus, differences in at lg 
and fig result in widely varying characteristics even 
for flames having nearly the same Lf. 


Turbulent Flames 

In turbulent non-premixed jet flames, D is not 
fixed but rather is nearly proportional to u'L { . Be- 
cause ti' - U 0 and L { - d 0 , L f - UJ*/u'Lt - d 0 for 
round jets; thus, L f is independent of U n . This pre- 
diction is supported by classical lg experiments [.53] 
as well as recent fig experiments [54] (Fig. 11). Note 
that Lf(jig)/Lf(lg) is practically constant even be- 
yond the transition to turbulent conditions (high U (> 
and Re). Note also that the maximum Re at which 
flame exists (the “blow-off* limit) is different at lg 
and fig. This is surprising because blow-off condi- 
tions are typically controlled by behavior near the 
flame base [55], where buoyancy effects are often 
considered insignificant. This suggests that blow-off 
is partially affected by convection induced by the 
buoyant plume far above the jet exit, even at very 
high U Q . One would intuitively conclude that the lg 
flames should blow-off at lower U n because buoyant 
flow would induce higher "effective** U w which is 
consistent with Fig. 11. This shows that buoyancy 
effects are quite ubiquitous even under conditions 
commonly thought to be unaffected by buoyancy. 
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Fig. 12. Direct photographs of sooting n-C 4 H 10 non- 
premixed gas-jet flames at lg (left) and^g (right) at Re *■ 
42, jet diameter 10 mm, showing evidence of thermopho- 
resis-indueed agglomeration at fi g. Photographs courtesy 
of Professor O. Fujita. 

Soot Formation Processes 

Compared with lg flames fig gas-jet flames have 
much greater tendencies to emit soot [49,56], indi- 
cating that increases in f Jet (thus, greater time for 
soot formation) plus broader regions in which com- 
position and temperature are favorable for soot for- 
mation [56,57] outweigh lower temperatures at fig, 
which decreases soot formation [58] . Recent quan- 
titative measurements [56] show peak soot-volume 
fractions about twice as high at fig than lg for 50% 
C2H2/50% N^-air flames. 

Surprisingly, fig gas-jet flames exhibit “smoke 
points,” corresponding to critical U 0 below which 
soot is consumed within the flame, and above which 
soot is emitted from the flame [57]. Smoke points 
are expected for buoyant round-jet flames because 
fj et ~ U^ 2 \ thus, increasing U Q increases the time 
available for soot formation. But for nonbuoyant 
flames, ~~ L™, suggesting no smoke point should 
exist. Fully elliptic numerical computations [59] 
show that, for some circumstances, fj et does increase 
monotonicallv with U ot which coula explain smoke 
points for nonbuoyant flames. This behavior was 
suggested to result from axial diffusion effects [59], 
but in this case, #j et should asymptote to constant 
values for large Lp, where axial diffusion is negligible. 
Thus, simple explanations of fig smoke points remain 
elusive. Residence-time considerations alone may be 
misleading; soot precursor temperature-composi- 
tion-time history effects are apparently also impor- 
tant. 

With weak convection, thermophoretic forces, 
which move particles toward lower temperatures, 
are an important effect. If convection and tempera- 
ture gradients are in the same direction, the convec- 
tive and thermophoretic forces may balance at some 
location. This leads to soot accumulation inside the 
flame front at fig (Fig. 12) [60], where an annulus of 


soot accumulation is convected by axial flow- through 
the flame tip. Downstream, for reasons not yet ex- 
plained, the soot annulus fragments, creating crown- 
like structures. These effects are only observable at 
fig where convection velocities are comparable to 
thermophoresis velocities, (** 5 mni/s for the con- 
ditions of Fig. 12 [60]). 

Condensed-Phase Combustion 

Droplet Combustion 

The first microgravity combustion experiments 
were isolated fuel-droplet tests conducted by Ku- 
magai and Isoda [61]. At lg. experimental measure- 
ments of droplet burning rates are compromised by 
buoyant convection, which destroys the spherical 
symmetry and inducing additional heat and mass 
transport, which alters burning rates and compli- 
cates modeling. Classical theory [62,63] predicts that 
burning rates for spherically symmetric buoyancy- 
free droplet burning are given by 

d%, ~ dl = ft; 

K - (8A//><jC p ) ln(l + B) (8) 

In experimental studies, the droplet diameter 
could be fixed by forcing fuel tnrough a porous 
sphere at rates that balances evaporation, leading to 
a steady mass burning rate (m) = {it/A)p^d ( \K, but 
most experiments employ fuel droplets where d ( \ de- 
creases with time. Somewhat surprisingly, many 
fuel-droplet results follow equation 8 well despite 
unsteadiness, heat losses, soot formation, and water 
absorption effects discussed below. 

As with flame balls, steady solutions exist for 
spherical droplet flames even in infinite domains, 
with the flame front located at 

d[ — d d ln(l 4- B)/ln(l + /) 

Whereas flame balls are convection-free, droplet 
flames exhibit Stefan flow due to fuel vaporization 
at the droplet surface. Mass conservation dictates 
that the Stefan velocity decays as 1/r 2 , causing tem- 
perature and concentration profiles to vary with ra- 
dius in proportion to (1 4- rather than 1/r 

as in flame balls, although these profiles, when nor- 
malized by flame radius, are indistinguishable at 
large r. Unlike flame balls, heat losses are not re- 
quired for stable droplet flames because the flame 
location cannot move away from the stoichiometric 
contour. 

The characteristic time scale for droplet combus- 
tion is fdrop df /a. This leads to two extinction lim- 
its [64,65], one for small df where f ( j rnp < t c \ wm and 
thus fuel and oxidant cannot react before interdif- 
fusing, and one for large d( where t^rop > Cad and 
thus the temperature decrease from radiative loss 
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Fig. 13. Effect of initial droplet diameter on quasi- 
steadv burning rate ( K ) for heptane droplets burning in air 
and an CL-He atmosphere a t jug, showing that K decreases 
with increasing apparently because of increased ac- 
cumulation of soot and gas-phase radiant species for larger 

4b m. 



Fig. 14. Direct photographs of heptane droplets burning 
in air at jig showing spherically symmetric combustion 
(left) and a soot “tail” formed by weak convection effects 
(right) [68], Photographs courtesy of Professor T. Avedi- 
sian. 


reduces the reaction rate sufficiently to cause ex- 
tinction. The former occurs when d( < (af c hem) 1/2 
and the latter when d( > Recent space ex- 

periments have reported radiative extinction of large 
droplets in air [66] and CL-He [67] atmospheres. 
Even when radiation does not cause extinguishment, 
it causes a decrease in K , especially when soot for- 
mation is significant [68,69]. 

In experiments, K is not constant as the quasi- 
steady model (equation 8) predicts, but rather de- 
creases with increasing d do (Fig. 13) [70]. This sug- 
gests that nonsteady effects, specifically the diffusion 
of thermal energy and radiant combustion products 
into the far field, may be significant; this would cause 
changing radiative loss over time. Analogous behav- 
ior occurs in flame balls [41], where radiative loss 
requires ~ 100 s to reach steady state even in flame 


balls much smaller than typical droplet flames. An- 
other indication of unsteadiness in droplet flames is 
that constant df/d { j values are not generally achieved, 
especially for large droplets at /ig [69, 70] /in contrast 
to the predictions of the quasi-steady theory. Un- 
steadiness effects in droplet flames were recentlv an- 
alyzed by King [71]. The effect of d tln on K has' also 
been proposed [68] to result from soot accumula- 
tion, whicn is more significant for larger </,i„ and acts 
to decrease net heat release and increase radiation. 

As with non-premixed gas-jet flames, soot particles 
in jug droplet flames exhibit thermophoresis effects 
[72,73], leading to soot agglomeration between the 
droplet and flame front (Fig. 14). The agglomerates 
may break apart suddenly, leading to multiple burn- 
ing fragments. Because velocity and temperature 
gradients are readily modeled in spherical droplet 
flames, such experimental observations enable as- 
sessment of thermophoresis effects on soot particle 
transport. 

Another complicating factor arises in fuels that are 
miscible in water. The fuel may absorb water vapor 
from the combustion products, causing significant 
departure from equation 8 and reducing flammabil- 
ity. Such behavior is found [74,75] in methanol 
flames, where extinction diameters depend substan- 
tially on d do due to water absorption during com- 
bustion. 

Finally, in recent space experiments, flame oscil- 
lations with amplitudes comparable to the mean 
flame diameter nave been observed [76]. The oscil- 
lation amplitude grew with time with extinction typ- 
ically occurring after eight cycles. The oscillation fre- 
quency was approximately 1 Hz. Cheatham and 
Matalon [65] predicted oscillations of roughly the 
correct frequency in droplet flames of pure fuels un- 
der near-extinction conditions when the reactant 
Lewis numbers are sufficiently high. To date, oscil- 
lations have been reported for methanol/dodecanol 
(80/20 mass fraction) bicomponent droplets burning 
in air; it is uncertain whether oscillations also occur 
in pure fuels. Moreover, jig experiments in 0 2 -He 
atmospheres [67], with much higher Le, did not ex- 
hibit oscillations. One possible explanation is that 
droplet support fibers were used in Ref. [76] but not 
Ref. [67]. The fiber could increase conductive and 
radiative losses, which encourage oscillations [65]. 
Another explanation is that, because oscillations oc- 
cur only near extinction conditions, depletion of ox- 
ygen (leading to extinction) was much more signifi- 
cant in the smaller combustion chamber used in the 
methanol/dodecanol experiments. Such oscillatory 
instabilities are discussed further in the following 
section. 

Candle Flames 

An excellent example of the differences between 
lg and jig flames is seen in perhaps the most com- 
mon and familiar of all combustion processes — can- 
dle flames. At lg, candle flames are supported by air 
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candle flames were noted [65 j. Alternatively, Buck- 
master [79] showed that the flame “edge" separating 
burning and nonbuming regions of non-premixed 
flames exhibit oscillatory behavior at (for quasi-sta- 
tionary edges) Le > l 4- S//?(l - { TJT f .)) 2 when 

Le for the other reactant is unity’. Le ** 1 for 0 2 in 
N 2 , but Le for fuel vapors is probably closer to 2; 
thus, edge-flame instabilities could explain the ob- 
served oscillations. While neither instability mecha- 
nism has been definitively linked to the candle-flame 
experiments, both predict greater propensities for 
oscillation with greater heat losses, which is consis- 
tent with the observation that oscillations occur near 
extinction. 


Fic. 15. Direct photographs of candle Hames at lg and 
fig, showing impact of buoyant flow on flame shape [77], 


entrained via buoyant flow, which generates a self- 
sustaining flame and flow configuration. Obviously, 
this mechanism cannot apply at fig. The spherical 
diffusion equation admits steady solutions for flame 
balls and droplet flames without forced convection. 
An interesting question is whether candle flames, 
which are not strictly spherical, can behave similarly. 
Space experiments [77,78] indicate that candle 
flames can be steady for >45 min, with flame shapes 
typically hemispherical (Fig. 15). Eventually, the 
flames always extinguished, whereas the spherical 
flame model predicts that the flame would bum in- 
definitely. Protective screens used in the experi- 
ments may have limited the 0 2 supply, eventually 
allowing sufficient 0 2 depletion to cause extinguish- 
ment. 

Before extinction, the candle-flame edge fre- 
quently advanced and retreated periodically. The os- 
cillation amplitude increased over time, and on one 
retreating cycle, the entire flame extinguished. With 
larger wicks (thus, larger flame diameters), oscilla- 
tions started spontaneously, whereas with smaller 
wicks, oscillations occurred only when solid objects 
were placed near the flame. Only a few cycles before 
extinction were observed in the Spacelab experi- 
ments [77], whereas hundreds of cycles were ob- 
served in the Mir experiments [78]. This is probably 
because the protective screen was much more per- 
meable in the Mir experiments, thereby decreasing 
the 0 2 depletion rate and maintaining the flame at 
near-extinction conditions much longer. 

At least two possible explanations for these oscil- 
lations have been advanced Cheatham and Matalon 
[65] showed that, near extinction, oscillatory insta- 
bilities occur in spherically symmetric droplet flames 
with radiative loss at sufficiently high Lewis num- 
bers. Their predicted oscillation frequencies (0.7- 
1.4 Hz) are comparable to experimental observa- 
tions, however, the differences between spherically 
symmetric droplet flames and roughly hemispherical 


Flame Spread over Solid Fuel Beds 

Flame spread over solid fuel beds is typically clas- 
sified as opposed- flow, where convection opposes 
flame propagation, or concurrent -flow. Because up- 
ward buoyant flow opposes flame spread, lg down- 
ward flame spread is opposed-flow, whereas upward 
flame spread is concurrent- flow. At fig without 
forced flow, flame spread is always opposed-flow' be- 
cause the flame spreads toward the fresh oxidant 
with a self-induced velocity equal to the spread rate 
(Sf). At lg, self-induced convection is negligible be- 
cause buoyancy-induced flow's are typically (ga g ) 1/3 
~ 0.10 m/s Sf. Very few concurrent- flow flame 
spread have been conducted at fig [80], conse- 
quendy, this section focuses on opposed-flow spread. 

Sf is estimated by equating tne conductive heat 
flux to the fuel bed ( = X(SW)(Tf - T v )/S . where 
S = a IU is the thermal transport zone thickness and 
U is the opposed-flow velocity (forced, buoyant and/ 
or self-induced), to rate of fuel bed enthalpy increase 
(= P,C Fs t s (T v — TJWSf). Assuming mixing-limited 
reaction (infinite-rate chemistry), for thermally thin 
fuels, where heat conduction through the solid is 
negligible, Sf is predicted to be [81,82] 


n 7> - 7\ 

4 p s C p . s r s T v - T x 


(9) 


Note that Sf is independent of U and P. For ther- 
mally thick fuels, where heat conduction through the 
solid fuel dominates, r s is the thermal penetration 
depth into the solid, estimated by equating the con- 
ductive heat flax to the fuel bed to the heat flux 
through the fuel (= A sy (<JW)((T v — T x )/r s ), where 
the subscript y refers to the direction normal to the 
fuel surface). This leads to the exact solution [81] 


= 


K(T V - TJ 
MTf - TJ/S 

=> S r - U 


K(T V - TJ a 
A(T r - 7V) S f 
XpCp ( T( — T v 


^•sACp.s \TV T z 


( 10 ) 


Note that, unlike the thin-fuel case, for thick fuels 
S r ~ U l P l . 
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Fig. 16 Minimum mole percent CL in N 2 supporting 
flame spread over a thin solid fuel bed, as a function of the 
opposed flow velocity ((./) [83], showing dual-limit behav- 
ior, that is, residence-time limited extinction at high U and 
radiative loss extinction at low U. 


Dual-limit extinction behavior is observed in Mg 
flame-spread experiments (Fig. 16) [83]. The time 
for thermal energy to diffuse across the convection- 
diffusion zone (f.iift ) is S/U = a/U thus, high-U 
extinction occurs when < f chem or U > 
(a/f c li ein ) 1/2 and radiative extinction occurs when 
> t rj( j or U < (aJt r . d j) 1/2 . (Surface radiative loss may 
also be important, particularly at moderate and 
higher U [84].) Interestingly, the minimum 0 2 con- 
centration supporting combustion (^ 02 j im ), and thus 
the greatest hazard, corresponds to 1/ 88 0.1 m/s, 
which is lower than buoyant convection at lg, and 
might correspond to ventilation drafts in manned 
spacecraft. 

A radiative loss parameter can be defined a s H = 


^litf ^rad - o^rkI- Because H — U S f is lower at 
where U is lower and thus H is higher. Experi- 
ments [8.5] (Fig. 17) show that an imposed forced 
flow at //g increases S, because U (sum of forced flow 
and S f ) increases and thus H decreases, whereas at 
higher U (whether buoyant or forced), S f decreases 
as the high-t/ limit is approached. For 21% 0 2 or 
lower, the infinite-rate chemistry prediction of equa- 
tion 9, Sr — t/ 0 , is never achieved. Only at 30% 0-> 
is T f high enough that this condition is achieved. 

Because a 1 and f ra(1 - F\ H - P >, thus for 
thin fuels, Sf should increase with P toward the ideal 
(adiabatic) value (equation 9). This is confirmed bv 
quiescent thin-fuel space experiments [S6,8?f, 
which show S f increasing from 3.2 to 5.9 mm/s as P 
increases from 1.0 to 2.0 atm with fixed 0 2 mole 
fraction (0.50). For these conditions, H decreases 
from 24 to 3.5, thus even at the highest P , radiative 
effects are probably still important. This is consistent 
with computations [86] that predict Sf = 12 mm/s 
(almost independent of P) for adiabatic conditions 
for this fuel/atmosphere combination. 

Neither N 2 nor O? emit thermal radiation: thus, 
for flames in 0 2 -N 2 atmospheres, only H-.0 and 
C0 2 combustion products radiate significantly. For 
these cases, ripically a? 1 1.2 m » S and thus, 

radiative transfer is optically thin (negligible reab- 
sorption). When a^S 5: 1, reabsorption effects can- 
not be neglected. With reabsorption, some radiation 
is not lost and may augment conduction to increase 
Sf above that without radiation. This behavior is seen 
experimentally [88] using strongly emitting/absorb- 
ing C0 2 and SF 6 diluents (Fig. 18), where S r is 
higher and/ 02 )im is lower at/ig than lg, whereas the 
opposite (conventional) behavior is found in nonra- 
diant diluents (not shown). These data indicate that, 
for nonradiating diluents. Mg is less hazardous be- 
cause *o 2 .iim is higher at Mg than at lg (0.21 vs. 0.16 
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Fig. 17. Flame-spread rate over a 
thin solid fuel bed as a function of 
the opposed-flow velocity (L r ), for 
three values of the mole percent 0 2 
[85 1, showing dual-limit response, 
Note that the infinite-rate kinetics 
prediction [81,82] that the spread 
rate is independent of U is only sat- 
isfied at 0 2 mole fractions higher 
than that in air. 
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FiC. 18. Flame-spread rate over a thin solid fuel bed at 
lg (downward spread) and fig as a function of O? mole 
fraction [88], showing atypical behavior in O 2 -CO 2 and 
Oi-SF fi atmospheres where the spread rates are higher and 
the minimum 0 2 model fraction supporting combustion is 
lower at fig. 




Fig. 19. Fingering patterns observed in smoldering 
flame spread over a thin paper fuel sample in fig (91). 
Flaming combustion was inhibited by soaking the fuel sam- 
ple in potassium acetate. An imposed convective velocity 
of 0.05 m/s flows from right to left. Grid pattern scale is 10 
mm by 10 mm. Photograph courtesy of Dr. S. Olson. 
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for He), whereas for radiant diluents, fig is more 
hazardous (0.21 vs. 0.24 for C0 2 ). This is particularly 
significant considering that C0 2 -based fire suppres- 
sion systems will be used on the International Space 
Station. To date, flame spread calculations have em- 
ployed optically thin radiation models with constant 
fl P [86,87] or variable depending on local tempera- 
ture and composition [80] and thus cannot assess 
reabsorption effects. 

These discussions pertain to thin fuels, for which 
steady jug spread is possible because theoretically Sy 
— L m . For thick fuels, Sy ^ U l , thus Sy is indeter- 
minate for quiescent /ig conditions ((/ = Sy). When 


unsteady solid-phase conduction is considered. r v — 
(a N t) l/2 , which results in Sf — t“ 12 [89). Conse- 
quently, all fuel beds at quiescent fig conditions 
eventually become thermally thin (penetration 
depth greater than the bed thickness) unless the ra- 
diative effects, discussed later, are considered. Of 
course, flames may extinguish due to large 6 (thus 
large radiative loss) before reaching steady-state, 
thermally thin conditions. 

A difficulty in comparing space experiments to 
two-dimensional model predictions is that the fuel 
bed width W (30 mm for thin fuels [86,87] and 6.2 
mm for thick fuels [89]) is smaller than the thermal 
transport zone thickness (S). Consequently, these ex- 
periments can hardly be considered two-dimen- 
sional. Both lateral heat loss, which retards spread, 
and lateral 0 2 influx, which enhances spread, are 
probably important, thus their effects may partially 
cancel. Some authors [86] suggest that radiative 
losses decrease S to values much smaller than aJU 
= a/Sy, but the oxygen transport zone thickness 
(Sq 2 ) is still D 0 /U, because no analog to radiative 
loss exists for 0 2 transport. Because Le = a/D 0 , ** 
l.*o. is nearly the same as S in adiabatic flames. 
Thus, the/ig flames mentioned above have probably 
benefited substantially from lateral 0 2 influx, espe- 
cially for lower pressures and 0 2 mole fractions, 
where S/\V ~ a/SyW is largest. In fact, Sy might be 
higher at smaller W due to lateral 0 2 influx. Space 
experiments using cylindrical fuel rods are planned 
[90] to examine truly two-dimensional spread. 

Recently, a surprising observation of fingering 
fronts was found in space experiments using paper 
samples treated to inhibit flaming combustion but 
allow smoldering propagation (Fig. 19) [91]. Finger- 
ing was observed at fig when U < 50 m m/s, whereas 
smooth fronts were observed at lg for all U. This 
was proposed [91] to result from limited 0 2 mass 
transport at fig with low U, which caused the 0 2 
consumption regions to become localized spots in- 
stead of continuous fronts. This proposition does not 
explain why heat conduction does not smooth out 
potential fingers as it does (for example) in pre- 
mixed-gas flames with Le > 1. The following alter- 
native explanation is proposed here. Gas-phase heat 
transport occurs on the length scale 5 ~~ a/U, and 
solid-phase transport occurs on the scale S s ~ otju^ 
where u s is the smolder front velocity and a, the 
solid thermal diffusivity. Oxygen transport occurs 
only through the gas phase on the scale D 0 JU 
a/JJ ~ S. Radiative loss can suppress heat transport 
through the gas, but no corresponding effect on 0 2 
transport can occur. Thus, at low U, the effective Le 
is a s /D Ql C 1. At higher U or at lg, S is smaller, gas- 
phase heat transport dominates, and radiative effects 
are weaker; thus, the effective Le is a/D 0 , ** 1. 
These assertions are consistent with estimates [91] 
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Time (s) 

Fig. 20. Measured (thick lines) and computed (thin 
tines) Hame position versus time for flame spread over a 1- 
butanol pool 20-mm wide and 25-mm deep [97]. The 10. 
20, and 30 notations refer to the opposed-flow velocities 
([/) in cm/s. Both computed results are for pg conditions, 
V = 30 cm/s, either with or without hot gas expansion 
(HGE). The comparison of predicted and measured results 
suggests a very strong influence of expansion, which is 
much less effective in the experiment because of the re- 
laxation of expansion in the transverse dimension, a factor 
not captured within the two-dimensional model. 

of the relative importance of gas-phase and solid- 
phase transport. Both premixed [27,39] and non- 
re mixed [92,93] flames with effective Le < 1 ex- 
ibit diffusive-thermal instabilities that cause 
fingering patterns, whereas for Le 2: 1, the fronts 
are stable. This explanation is also consistent with Ig 
experiments [94] on horizontal fuel beds burning in 
oxidant channels of adjustable vertical height At 
small heights or low U, fingering similar to Fig. 19 
was observed. In this case, conductive loss to the 
channel ceiling causes suppression of gas-phase heat 
transfer. Apparently, in both cases, the key factor is 
suppression of gas-phase heat transfer while allow- 
ing solid-phase heat transfer, which reduces the ef- 
fective Le (though this factor was not mentioned in 
Ref. [91] or Ref. [94].) 

Flame Spread over Liquid Fuel Pools 

Flame spread over liquid fuels encompasses prac- 
tically all solid-fuel flame spread phenomena dis- 
cussed above, plus liquid-phase flow effects. Typi- 
cally, T v — T x is smaller for liquid than solid fuels, 
thus Sf is higher. Also, if T v - T* is small, some fuel 
prevaporization occurs even at T — T x , thus partially 
premixed gas-phase combustion phenomena may 
occur Because of the fuel surface temperature gra- 
dient upstream of the flame, surface tension gradi- 
ents are produced that cause the surface layer to 


move upstream (away from the flame), which in- 
creases Sf. At lg, this heated liquid layer must lie 
near the surface, whereas at /ig, no limitation exists. 

Experiments at lg, suminarizedfin [95], show that 
at low fuel temperatures, the average Sfis small (typ- 
ically 10 mm/s) and the spread alternates between a 
fast “jump ' velocity' and a slow “crauT* velocity. At 
higher fuel temperatures, Sf is faster and steady. For 
the conditions exhibiting pulsating spread at lg, fig 
flame spread cannot be maintained, whereas for the 
conditions exhibiting uniform spread at lg, steady 
spread is also exhibited at pg [96]. Pulsating spread 
has never been observed at pg. No definitive expla- 
nation for these obse nations has been advanced. 
For the conditions exhibiting pulsating spread at lg 
and no spread at pg, flame spread is still different at 
lg and pg when forced flows comparable to buoy- 
ancy-induced convection (U — 0.30 m/s ) are im- 
posed (Fig. 20) [97]. Specifically, lg spread is almost 
unaffected by the imposed flow, but pg spread is 
steady with Sf being lower (**15 mm/s) than either 
the lg jump velocity- (^ 100 mm/s) or crawl velocity 
(^*22 mm/s). 

Detailed numerical modeling [98] predicts pulsat- 
ing spread at pg for the conditions of Fig. 20 and 
values of Sf much closer to the measured lg Sf. Re- 
markably, if thermal expansion is artificially sup- 
pressed, good agreement between the model and pg 
experiments is found. It is proposed [97] that this 
agreement results from three-dimensional effects; 
specifically, in the experiment, flow induced bv ther- 
mal expansion is relaxed in the lateral dimension, 
whereas the two-dimensional model does not permit 
this. That three-dimensional effects might dominate 
is surprising considering that, for this flame, S/\V ~ 
a/UW *=* 0.02, thus, 3 C W. Also, this hypothesis 
does not explain why pulsating flame spread is ob- 
served at lg but not pg. Consequently, in the case 
of liquid-fuel flame spread, pg experiments have 
identified limitations in our current understanding 
of combustion processes at lg. 

Recommendations for Future Studies 

Re absorption Effects 

The pg studies described here suggest new unre- 
solved issues and opportunities for further improve- 
ments in understanding. Perhaps the most important 
are the effects of reabsorption of emitted radiation, 
including both reabsorption by the emitting gas and, 
in two-phase combustion, absorption by the con- 
densed phase. All radiative effects discussed above 
are critically dependent on the degree of reabsorp- 
tion. To study reabsorption effects requires radia- 
tively active diluents (C0 2 , SF 6 ), high pressures and/ 
or large systems. Ail of these conditions lead to 
higher Gr d at lg and thus turbulent flow. Hence, pg 
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experiments enable study of reabsorption effects 
without the additiomil complications due to turbu- 
lence. 

Reabsorption effects are important not only to fig 
studies but also to combustion at high pressures and 
in large combustors. For example, at 40 atm, typical 
of premixed-charge internal combustion engines, a P 
** ^ ni 1 ' thus, a p 1 — 0.045 in, for stoichiometric 
combustion products. Because this length scale is 
comparable to cylinder radii, reabsorption effects 
within the gas cannot readily be neglected. Simple 
estimates [99] indicate radiative loss may influence 
flame quenching by turbulence in lean mixtures. 
Similarly, reabsorption cannot be neglected in at- 
mospheric-pressure furnaces larger thanflp 1 ** 2.2 
m. Moreover, many combustion devices employ ex- 
haust-gas or flue-gas recirculation; for such devices, 
the unbumed mixtures contain significant amounts 
of absorbing C0 2 and H 2 0. 

Although reabsorption of emitted radiation could 
affect practically all types of flames reviewed here, 
to date, reabsorption effects have been studied only 
for propagating pre mixed-gas flames [24,26], flame 
balls [45], and flame spread over thermally thin fuels 
[88]. All have shown substantial differences from op- 
tically thin behavior. Two examples of effects ex- 
pected for other flames are given below. 

Reabsorption effects could be substantially more 
important for droplet combustion than for flame 
spread over solid fuels because for droplets, the Ste- 
fan flow severely limits heat conduction to the drop- 
let surface. This is why heat release (B) affects burn- 
ing rates ( K ) only weakly (logarithmically) (equation 
8). Radiative transfer is unaffected by the Stefan 
flow. Equation 8 is readily extended to include sur- 
face radiative flux ( q r ): 


O — ln( 1 + 


B 


1 - R/Q 


__ q4dPp 

* 


. KpdPp 
8/1 
(ID 


a fd t ^2)}o(T A Ti) [102]. Using volume-averaged 

properties T = 1000 K andrt P( = 20 m R = 
18 and thus = 8.0, indicating radiation 

dominates heat transport. As discussed later, at high 
pressures radiative effects may prevail even in CK~ 
N 2 atmospheres. 

Flame spread over thermally thick fuel beds in 
quiescent atmospheres at fig is typically considered 
inherently unsteady [89], however radiative transfer 
to the fuel bed could enable steady spread. If the 
flame is modeled by an isothermal volume with di- 
mensions 6 X S X W, radiation induces a radiative 
flux A6-W = A {a 2 /Sf)\V that augments the con- 
ductive flux k(SW)(Tf - T k )/5 . Equating this total 
heat transfer to p s C Ps x^(T v - TJ WS, yields 

Sf = f — ±4 F 

W<ty<CpMT v - rj - ;jT f - rj- 

( 12 ) 

which vanishes without gas radiation (A = 0). Thick- 
fuel space experiments in 0 2 -C0 2 or 0 2 -SF 6 at- 
mospheres could be employed to check for steady 
spread and test the accuracy of equation 12. 

Whereas optically thin radiation modeling is rea- 
sonably straightforward, modeling of spectrally de- 
pendent emission and absorption is challenging be- 
cause local fluxes depend on the entire radiation 
field, not just local scalar properties and gradients. 
Some relevant computations have emploved gray- 
gas models [74], but recent studies [26,103] show 
that these methods are probably inaccurate because 
of the wide variation in spectral absorption coeffi- 
cient with temperature, species, and wavelength. 
Comparisons of various radiative treatments for 
small one-dimensional non-premixed flames have 
been made [103], Comparisons for larger, multidi- 
mensional systems would be valuable. Moreover, re- 
cent studies of/ig soot formation [52,56,57] may en- 
able improved modeling of soot radiation at fig. 


Figure 21 shows the predictions of equation 11. (Al- 
though equation 11 has apparently not been pre- 
sented previously, numerical studies [100,101] have 
shown qualitatively similar predictions. Moreover, 
these studies show that typical radiative absorption 
lengths for liquid fuels at relevant wavelengths are 
on the order of 1 mm, thus large droplets could ab- 
sorb most incident radiation.) For spherical shells of 
radiant combustion products having thickness 3 
df, q r Adi 4, then for typical values B = 8.5, A = 
2 x 10 6 Wl m 3 , S = 10 mm, d d = 5 mm, C y = 
1400 J/kgK, A = 0.07 W/mK, and L v = 400 kf/kg, 
equation 11 predicts R = 0.63 and Q/Q Rm0 = 1.11; 
thus, moderate effects of radiative transfer are ex- 
pected. For droplets in radiatively active diluents 
such as C0 2 , the effect could be much stronger. Us- 
ing the PI approximation, for a sphere of unit 
emissivity in an infinite gray gas, q T = [4/(2 4* 3 


High-Pressure Combustion 

All practical combustion engines operate at pres- 
sures much higher than atmospheric. The impact of 
buoyancy for premixed flames scales as - 

(ga/Sl) 2 ** - P n ~ 4n t where n is the overall reaction 
order (S L — P n/2 ~ i ). Because, typically, n < 4/3 for 
weak mixtures [104], where buoyancy effects are 
most important, the impact of buoyancy increases 
with pressure. Also, as discussed earlier, radiation 
effects are more difficult to assess at higher pressure 
due to increased interference from buoyant trans- 
port. Nevertheless, few high-pressure fig combus- 
tion experiments have been performed. High-pres- 
sure droplet combustion experiments [105] revealed 
substantial but different increases in K with P at lg 
and fig. Radiative effects were not discussed but 
could have been important, because in equation 11, 
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Fig. 21. Predicted effect of radiative heat transport co- 
efficient (R) on droplet burning rate constant referenced 
to the value without radiative transport, showing impor- 
tance of absorption of radiation at the droplet surface on 
the resulting burning rate. B = 3 and B = 8.5 are char- 
acteristic of methanol and heptane, respectively, burning 
in air. 

the only pressure-dependent factor i s q r — A ~ P l , 
thus in Fig. 21, R — P. (For most flames, length 
scales decrease with increasing P, which would de- 
crease radiative effects; but for droplet flames, df 
depends onlv on stoichiometry [75,76].) Conse- 
quendy, further assessment of radiative effects in 
high-pressure droplet combustion and other types of 
flames appears warranted. 

Three-Dimensional Effects 

In earlier sections, effects of lateral heat and mass 
transport on flame spread were discussed. To assess 
three-dimensional effects, fig experiments with vary- 
ing fuel bed width (W) are needed. Complementary 
three-dimensional modeling using codes such as 
those developed by NIST [106], extended to include 
gas-phase radiation, would be instructive. An ap- 
roximate but much less expensive approach would 
e to incorporate volumetric terms 6MT(x , y) - Tf)/ 
W 2 and BpDjYJx, y) — Y, JAV 2 , where x and y are 
the coordinates parallel and perpendicular to the 
fuel bed, into the two-dimensional model to account 
for lateral heat losses and lateral diffusion of each 
species i. 

Another three-dimensional effect is found in the 
development of flame balls from ignition kernels. 
Currently, it is known that large flame balls are lin- 
early unstable to three-dimensional disturbances for 
weak loss (Fig. 5), but the transition from splitting 
flame balls to stable flames is not well understood 
nor can the number of flame balls produced from an 
ignition source be predicted. Modeling using three- 
dimensional premixed flame codes [107] is needed. 

Gas-Jet Flames 

Table 2 shows predicted scalings of flame lengths 
(L f ) and residence times (f jet ) for buoyant and non- 
buoyant round-jet and slot-jet flames. Despite nu- 
merous investigations of round-jet flames at fig, no 


fig slot-jet results are available to test those predic- 
tions. Currently, it is unknown whether slot-jet 
flames at fig would exhibit smoke points or whether 
tfiis information could be used to explain smoke 
points in fig round-jet flames. Because L f depends 
on g for buoyant slot-jet but not round-jet flames, L { 
should be quite different at lg and fig for slot-jets 
but not round-jets. Residual accelerations in aircraft 
fig experiments will be more problematic for slot-jet 
than round-jet flames because f jH — g- 1/2 for buoy- 
ant round-jet flames whereas f i( . t — g- for slot-jet 
flames. 

There has been little investigation of blow-off be- 
havior of laminar gas-jet flames at/ig. Dual-limit be- 
havior might occur for flames of fixed- mass flow rate 
but varying d ()y with short residence time extinction 
at small d () (t^ t ~~ drJUJ and radiative extinction at 
large d n (thus, large f jet ). Experiments should be con- 
ducted by diluting the fuel rather than increasing U (> 
to obtain blow-off without transition to turbulence. 
In this wav, dual-limit behavior has been observed 
at Ig [92] with short residence time and conductive 
loss (to the burner rim) extinction branches. 

Quasi- Steady Spherical Diffusion Flames 

As discussed previously, comparing predicted ra- 
diative extinction limits of droplet flames to experi- 
ments is problematic because quasi-steady condi- 
tions may not be obtained, because extinction occurs 
for sufficiently large droplet and flame diameters, 
but the droplet diameter decreases throughout its 
life. Numerical models can account for transient ef- 
fects, but the multidimensional ignition process is 
difficult to model quantitatively. Comparisons of 
droplet experiments and computations to corre- 
sponding results obtained using fuel-wetted porous 
spheres would be most interesting. The fuel should 
be forced through the porous sphere at slowly in- 
creasing rates until extinction occurs when df > 
(<tf r ad) i/2 . thus obtaining truly quasi-steady extinc- 
tion. Lon gfig durations would be required to estab- 
lish steady diffusion-dominated far-field tempera- 
ture and composition profiles (thus, steady radiative 
loss). Some preliminary results have been obtained 
in aircraft experiments [108] but were severely com- 
promised by the residual g levels. Candle flames are 
similar to wetted porous spheres, though without 
true spherical symmetry or any means to control or 
measure the instantaneous d^ 

A simpler related experiment employs porous 
spheres through which gaseous fuel is forced at pre- 
scribed m, resulting in a flame diameter df ** mCp/ 
27rxln( 1 + /). Some drop-tower experiments using 
this configuration have been reported [109,110], al- 
though steady-state conditions were not obtained 
due to short fig durations. Steady-state behavior ap- 
pears unlikely given the m (thus, air consumption 
rates) and chamber sizes employed to date; near-wall 
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oxygen depletion would be significant over the times rates for H + O, + M reactions p irticul irlv the 
reqmred to reach steacK-state ( 2 10 s [ 1 10|). Steady- Chaperon efficienc.es of vanous vHecies f5l 
state behavior might be obtained in drop-towers bv These reactions ire extremelv i m l . ^ . ' 

using smaller m and diluted fuel with enriched ox- li.Jt flames due to 1 C* .JIT 

andff^-^/a'' t0 " ltreaSe fa " d thus decrease d f branching and chain-inhibiting steps near limits [21 1 
,r " ,> / but are much less important away from limits. Fur- 

ther scrutiny of the rates for these reactions at in- 
termediate temperatures (1100-1400 K) would be 
welcomed. 


Cataljtic Combustion 

Catalytic combustion holds promise for reduced 
emissions and improv ed fuel efficiency in many com- 
bustion systems [111,112], Because catalysis occurs 
at surfaces, it is inherently multidimensional and/or 
unsteady, requiring reactant transport to the surface 
and heat and pioducts transport from the surface. 
Although boundary -layer approximations can be in- 
voked, probably the only truly one-dimensional 
steady' catalytic configuration is a spherical surface 

immersed in nonbuoyant quiescent premixed gas 

a “catalytic flame bail.” In this case, r f is fixed but 
the surface temperature T s and fuel concentration 
Y s are unknown. These are related through the en- 
ergy conservation (including surface radiation) and 
diffusion equations to obtain the surface reaction 
rate in moles per second (O): 


9( Y s , T s ) = p s D,r s Y x ( 1 
T - T x 


Y, 

t ml ~ u 


T«j - T, 


Y/YJ/M 

o e.r,(T* - T*J 


Mi + 


a,(T, 




(13) 


where e, is the surface emissivity and properties with 
the subscript s are evaluated at T = T,. By varying 
r„ Y„ pressure, and diluent gas, &(T, Y) can be in- 
ferred from equation 13 and the measured T s . Of 
course, conditions must be unfavorable for initiation 
of propagating flames or flame balls that stand off 
from the surface. 


Chemical Models 

An important contribution of pp combustion ex- 
periments has been improved understanding of ex- 
tinction processes, which are inherently related to 
finite- rate chemistry. To obtain closure between ex- 
periments and computations, accurate chemical 
models are needed. For lean premixed hydrocarbon- 
air flames, most models [21,26] predict higher S L 
and leaner flammability limits than experimental ob- 
servations [22,33,113], The discrepancy seems larger 
than experimental uncertainty or unaccounted heat 
losses could explain. In contrast, for H r air flame 
balls [45] and lg strained premixed H r air flames 
[114], these chemical models predict smaller balls, 
lower S L , and richer flammability limits than exper- 
iments. These chemical models predict S L in 
mixtures away from extinction limits very faithfully. 
The discrepancies result largely from differences in 


Conclusions 

Our understanding of combustion fundamentals 
has been broadened by fig experiments into regimes 
not previously explored. In particular, these experi- 
ments have helped integrate radiation into flame 
theory. Although flame radiation has long been rec- 
ognized as an important heat-transfer mechanism in 
large fires [115], its treatment has largely been ad 
hoc because of the difficulty of predicting soot for- 
mation. Also, large-scale fires at lg are inevitably tur- 
bulent, leading to complicated flame-flow interac- 
tions Small-scale fig flames are laminar, often 
soot-free, and have significant influences of radia- 
tion. As a result of radiation effects, both premised 
and non-premixed flames have exhibited dual-limit 
extinction behavior, with residence time-limited ex- 
tinction at high strain, or curvature and radiative 
loss-induced extinction at low strain or curvature. 
The high-strain limit is readily observed at lg; when 
forced flow is absent, buoyant flow causes this strain. 
For weak mixtures, these limits converge, but the 
convergence and the entire low-strain extinction 
branch can only be seen at fig. This dual-limit be- 
havior is observed for stretched and curved pre- 
mixed-gas flames, strained non-premixed flames, 
isolated fuel droplets, and flame spread over solid 
fuels. Besides radiative effects, fig studies have en- 
abled observation and clarification of numerous 
other phenomena, for example, thermophoresis ef- 
fects in soot formation, spherically symmetric drop- 
let burning, diffusion -con trolled premixed flames 
(flame balls), and flame instabilities in droplets and 
candle flames. Considering the rapid progress made 
recently, further advances are certain to occur. It is 
hoped that this report on the current state of un- 
derstanding will help motivate and inspire such ad- 


vances. 


a P 

A 

B 


Cs 

C F 


Nomenclature 

Planck mean absorption coefficient 
flame surface area 

transfer number (chemical enthalpy gen- 
eration/enthalpy needed for fuel evapo- 
ration) 

stoichiometric molar ratio of fuel to air 
constant-pressure heat capacity 
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df 

<k 
d, |„ 
d„ 

D 

E 


Cr, 


H 


K 

L f 

L { 

L v 

Le 


m 

M 

P 

r 

l 


R, 

Re 

Sf 

^Uiin 

Sc 

^chem 

fditr 

^drop 

^inv 

^jet 

had 

his 

T 

T* i 
u' 

U 

u u 


Uo 

w 

\v 

y 

Y 

a 

fi 


characteristic flow-length scale or tube 
diameter 

droplet flame diameter 
droplet diameter 
droplet initial diameter 
jet exit diameter ( round jets) or slot width 
(slot jets) 
mass diffusivity 

overall activation energy of the heat-re- 
lease reactions 

stoichiometric fuel to air mass ratio 
acceleration of gravity 
earth gravity 

Grashof number based on characteristic 
length scale (d) = gdVv 2 
heat transfer coefficient in a cylindrical 
tube = 16 i/d 2 

radiative loss parameter for flame spread 

“ fdin^rat! = 

droplet burning rate constant (equation 8) 
flame length for gas-jet flame 
turbulence integral scale 
latent heat of vaporization of liquid fuel 
Lewis number (a/D = thermal diffusivity/ 
reactant mass diffusivity) 
mass burning rate 
fuel molecular weight 
pressure 

radial coordinate 
flame radius 

scaled flame radius (eq 6); radiation pa- 
rameter (eq 11) 
gas constant 

jet Reynolds number = UjiJv 
flame spread rate over solid fuel bed 
premixed laminar burning velocity 
burning velocity at the flammability limit 
Schmidt number = v/D 
chemical time scale 
thermal diffusion time scale = a/C/ 2 
droplet flame time scale = dj/a 
inviscid buoyant transport time scale 
residence time of non-premixed jet flame 
radiative loss time scale 
viscous buoyant transport time scale 
temperature 

adiabatic flame temperature 
turbulence intensity 
convection velocity 

local axial velocity in counterflow configu- 
ration 

jet exit velocity 

gas-jet flame width 

solid fuel bed width 

axial coordinate 

fuel mass faction 

thermal diffusivity 

nondimensional activation energy = 

E/RJ f 


/l >2. Inn 

6 

7 


P 

a 

L 


minimum oxygen concentration support- 
ing combustion 
flame thickness 
gas specific heat ratio 
thermal conductivity 

radiative heat loss per unit volume = 
4 oa P (T* - V.) 
kinematic viscosity 
density 

fuel bed half-thickness (thin fuel) or ther- 
mal penetration depth (thick fuel) 
Stefan-Boltzmann constant 
flame stretch rate 


Subscripts 


d 

f 

s 

V 


droplet surface condition 
flame front condition 
solid-fuel or solid-surface condition 
solid- or liquid-fuel vaporization condition 
- ambient conditions 
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COMMENTS 


H. S. Mttkunda, Indian Institute of Science. Bangalore, 
India. In your presentation, you credited the computa- 
tional demonstration of the result on the absence of a flam- 
mability limit in adiabatic flames to Ciovangigli and 
Smooke [1] This result was first obtained in 1988 (2J. You 
might like to amend the citation suitably. 
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Author s Reply. In this earlier paper by the questioner 
and his collaborators [2], it was shown that the solutions of 
the unsteady planar one-dimensional adiabatic premixed 
flame equations do not predict anv flammability limit for 
lean CH^Oo-NN mixtures — the burning velocity decreases 
asymptotically to zero as the mixture strength is decreased. 
In Smooke and Ciovangigli (Ref. [1] of the comment), the 


steady version of these equations was solved for CH 4 -air 
and H 2 -air mixtures, and a similar conclusion w as reached. 
Both are important findings, because it is possible in prin- 
ciple that unsteady effects could suppress limits of steady 
flames or cause limits to occur for mixtures that are flam- 
mable as steady flames. Together thev show that one must 
identify loss mechanisms of the appropriate magnitude to 
explain flammability limits, which was one of the kev points 
of my presentation. 


J P Gore. Purdue University, USA. You seem to have 
drawn a conclusion based on vour time-scale analysis that 
radiation is not important at 1 g. For combustion under 
high-pressure conditions, for fires (large and small), and 
for till flames as far as high activation energy reaction steps 
are invoked, radiation can be potentially important. Could 
you please re-examine and clarify your conclusion regard- 
ing radiation not being important at 1 g? 

Authors Reply. Certainly radiation is an important 
mechanism of heat transfer in flames even at 1 g— large 


li 
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in< lu.striiil furnaces and boilers depend almost entirely on 
iMt i Lit i\ * - transter The same is true ol many lar^e accidental 
fi r ,. s .i-.-r.. in lutililin^s. mvr liquid iuel pools?, I In* message 
l intended hi convey was that the radiation transport time 
stale is longer than the buoyant transport time scale* at 
earth -navitv. and tlms radiation dot's not play an important 
mle on local Name structure and extinction conditions at 


earth *»ravitv compan*d to the role that but >\ ant convection 
plavs Even on the lartje scale, however, buoyancy is prob- 
able a more dominant transport mechanism at j ^ — while 
a building fire would be very dillerent with buoyancy but 
no radiation, imagine bow dilleieut it would lie with radi- 
ation but no buoyancy! 



